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It can be computationally advantageous to perform computer simulations in a Lorentz boosted 
frame for a certain class of systems. However, even if the computer model relies on a covariant 
set of equations, it has been pointed out that algorithmic difficulties related to discretization errors 
may have to be overcome in order to take full advantage of the potential speedup. We summarize 
the findings, the difficulties and their solutions, and show that the technique enables simulations 
important to several areas of accelerator physics that are otherwise problematic, including self- 
consistent modeling in three-dimensions of laser wakefield accelerator stages at energies of 10 GeV 
and above. 



I. INTRODUCTION 

In 1], we have shown that the ratio of longest to 
shortest space and time scales of a system of two or 
more components crossing at relativistic velocities is 
not invariant under a Lorentz transformation. This 
implies the existence of an "optimum" frame of refer- 
ence minimizing a measure of the ratio of space and 
time scales. Since the number of computer operations 
(e.g., time steps), for simulations based on formula- 
tions from first principles, is proportional to the ra- 
tio of the longest to shortest time scale of interest, 
it follows that such simulations will eventually have 
different computer runtimes, yet equivalent accuracy, 
depending solely upon the choice of frame of refer- 
ence. The scaling of theoretical speedup was derived 
for a generic case of two crossing identical rigid parti- 
cle beams, and for three particle acceleration related 
problems: particle beams interacting with electron 
clouds Q, free electron lasers (FEL) Q, and laser- 
plasma accelerators (LWFA) 0] ■ For all the cases con- 
sidered, it was found that the ratio of space and time 
scales varied as 7 2 for a range of 7, the relativistic fac- 
tor of the frame of reference relative to the optimum 
frame. For systems involving phenomena (e.g., parti- 
cle beams, plasma waves, laser light in plasmas) prop- 
agating at large 7, demonstrated speedup of simula- 
tions being performed in an optimum boosted frame 
can reach several orders of magnitude, as compared 
to the same simulation being performed in the labo- 
ratory frame. 

We summarize the difficulties and limitations of the 
method, the solutions that were developed, and the 
simulations that we have performed to date. We show 
that the technique enables simulations important to 
several areas of accelerator physics that are otherwise 
problematic. For the first time, it allows for direct 
self-consistent simulations of laser wakefield accelera- 
tor stages at 10 GeV and beyond using current super- 
computers in a few hours, while the same calculations 
in the laboratory frame would take years using simi- 



lar resources and are thus impractical. It also allows 
simulations of electron cloud effects in high energy 
physics accelerators (modeled so far with codes based 
on quasistatic approximations) using more standard 
Particle-In-Cell methods. This renders these types 
of simulations accessible to a wider range of exist- 
ing computer codes, alleviates the added complication 
due to pipelining when parallelizing a quasistatic code, 
and removes the approximations of the quasistatic 
method which may not be applicable in some situa- 
tions. For free electron lasers, the new technique offers 
the possibility of calculating self-consistently config- 
urations that are not accessible with standard FEL 
codes due to the limitations of the approximations 
that they are based on. Finally, the method may offer 
a unique way of calculating self-consistently, and in 
three-dimensions, coherent synchrotron radiation ef- 
fects which are of great importance in several current 
and future accelerators. 



II. DIFFICULTIES 

Even if the fundamental electrodynamics and par- 
ticles equations are written in a covariant form, the 
numerical algorithms that are derived from them 
may not retain this property and special techniques 
have been developed to allow simulations in boosted 
frames. As an example, we considered in [j| an iso- 
lated beam propagating in the laboratory frame at 
relativistic velocity. When applying the effect of the 
beam field on itself using the Newton-Lorentz equa- 
tion of motion, the contribution from the radial elec- 
tric field is largely canceled by the contribution from 
the azimuthal magnetic field. However, we showed 
that the so-called 'Boris particle pusher' Q (which 
is widely used in PIC codes), does make an approxi- 
mation in the calculation of the Lorentz force which 
leads to an inexact cancellation of the electric com- 
ponent by the magnetic component. The magnitude 
of the error grows with the beam relativistic factor 
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and in practice, it is unacceptably large for simula- 
tions of ultra-relativistic charged beams, where the 
cancellation needs to be nearly complete. The issue 
was resolved by changing the form of the Lorentz force 
term in the Boris pusher, and solving analytically the 
resulting implicit system of equations (see [j| for de- 
tails). 

An additional practical complication of numerical 
simulation in a boosted frame is that inputs and out- 
puts are often specified (or desired) in the labora- 
tory frame. For example, in LWFA simulations, laser 
and plasma parameters have to be transformed from 
the laboratory to the new relativistic boosted frame, 
so that the electromagnetic waves will be Doppler- 
shifted, and the background plasma, with higher den- 
sity, is now drifting. In the PIC code Warp [§], the 
laser is injected at a plane that is fixed in the labo- 
ratory frame and drifting in the boosted frame. Like- 
wise, the initial phase-space distribution of a particle 
beam is generally known in the laboratory. For calcu- 
lations in boosted frames of large 7, deriving the initial 
beam conditions at a given time can be easy if the ini- 
tial conditions are simple (e.g., initial Gaussian beam 
in vacuum), or more difficult and/or computationally 
costly if injecting the beam in a particle accelerator for 
example, where its longitudinal extent in the boosted 
frame can cover several lattice periods. In order to cir- 
cumvent this difficulty, a procedure was implemented 
in Warp which injects the beam through a transverse 
plane that is fixed in the laboratory, but drifting in the 
boosted frame, similarly to the laser injection method. 
Due to long range space charge forces, it is still nec- 
essary to provide a reasonable estimate of the beam 
distribution near the injection plane; this is accom- 
plished by the use of "frozen" drifting macroparticles. 

After the relativistic PIC algorithm evolves the sys- 
tem in the boosted frame, the results must be trans- 
formed back to the laboratory frame. We have found 
it convenient in Warp to record quantities at a number 
of regularly spaced "stations" , immobile in the labora- 
tory frame, at a succession of discrete times, for both 
detailed time histories and laboratory time-averages. 
Since the space-time locations of the diagnostic grids 
in the laboratory frame generally do not coincide with 
the space-time positions of the macroparticles and 
grid nodes used for the calculation in a boosted frame, 
some interpolation is performed during the data gath- 
ering process. 

Finally, in simulations of laser-plasma acceleration 
stages (see below), we observed a short wavelength 
instability with a growth rate that rises with the ve- 
locity of the boosted frame and the inverse of the grid 
resolution, which we have controlled through the use 
of low dispersion electromagnetic solvers [7| and low- 
pass digital filtering. The details of the instability and 
its cures will be detailed in a future paper Q . 

Together with mitigation of numerical artifacts as 
just described, these techniques allow simulations us- 



ing boosted frames, with orders of magnitude speedup 
over the same simulations performed using a labora- 
tory frame, as shown below. Additional details of the 
input and output procedures can be found in [lOfl . 



III. EXAMPLES OF APPLICATION 
A. Laser wakefield acceleration 

Laser driven plasma waves offer orders of magni- 
tude increases in accelerating gradient over standard 
accelerating structures (which are limited by electri- 
cal breakdown), thus holding the promise of much 
shorter particle accelerators. Yet, computer model- 
ing of the wake formation and beam acceleration re- 
quires fully kinetic methods and large computational 
resources due to the wide range of space and time 
scales involved [Hj]. For example, modeling 10 GeV 
stages for the LOASIS (LBNL) BELLA proposal [Hi 
in one-dimension demanded as many as 5,000 pro- 
cessor hours on a NERSC supercomputer [15| . As 
discussed in Q, the range of scales can be greatly 
reduced if one adopts the common assumption that 
the backward-emitted radiation can be neglected, en- 
abling, for the first time, the full-PIC simulation of 
the next generation of laser systems. 

Warp simulations at plasma density n e = 10 19 cm~ 3 
were performed in 2-1/ 2D and 3D using reference 
frames moving anywhere between 7^ = 1 (laboratory 
frame) and 10. These simulations are scaled replicas 
of 10 GeV stages that would operate at actual densi- 
ties of 10 17 cm~ 3 [l(| [13] and allow short run times 
to permit effective benchmarking between the algo- 
rithms. Agreement within a few percent was observed 
on the beam peak energy and average energy between 
calculations in all frames, showing that the boosted 
frame simulations gain speed without sacrificing ac- 
curacy. A speedup of 100 was measured between the 
calculation in the frame at 7 = 10 and the calculation 
in the laboratory frame. 

The boosted frame model was then used to con- 
duct full scale simulations of 10 GeV stages at plasma 
densities of 10 17 cm -3 in 2-1/2D and 3D simulations. 
Simulations at 10 18 cm -3 were also conducted to es- 
tablish scaling. Relativistic factors of the boosted 
frame were 130 and 40 respectively, i.e. close to the 
relativistic factor associated with the wake velocity 
in the laboratory frame. The 3D run at full scale 
took almost 4 hours with 7/ = 130 using 512 cores 
on the cluster Lawrencium at LBNL. This provided 
direct simulation of next generation experiments and 
possible laser-plasma collider stages. Good agreement 
with the scaled energy gain was obtained. Such sim- 
ulations are impractical in the laboratory frame, with 
projected time of 15 years on the same resources us- 
ing the 27^ formula for the estimated speedup, and 
scaling from standard PIC runs. 
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B. Electron cloud driven instabilities 



Several existing and planned future particle accel- 
erators have limitations due to the electron cloud in- 
stability that may negatively impact the beam quality 
and in some cases even lead to severe beam loss. A cal- 
culation of electron cloud driven instability [2J for an 
ultra-relativistic beam was performed with the Warp 
code framework in (a) standard PIC mode using the 
new particle pusher in a Lorentz boosted frame; (b) in 
quasistatic mode [ll[ using linear maps to push beam 
particles into the accelerator lattice. The two runs 
were in good agreement and completed using similar 
computer resources and runtimes. The speedup fac- 
tor of the PIC boosted frame calculation compared 
to a PIC calculation in the laboratory frame was esti- 
mated at 500. For many calculations of electron cloud 
instability, the boosted frame approach may not re- 
solve any additional physics not included in the qua- 
sistatic approach. We note, however, that the qua- 
sistatic method requires significant special coding to 
take into account eventual longitudinal motion of elec- 
trons [llj], as wen as a special parallelization scheme 
[T3 | for parallelization along the axis of beam propa- 
gation, which are not standard to PIC codes. By con- 
trast, the boosted frame method includes naturally 
the longitudinal dynamics and requires more mod- 
est modifications to an existing standard PIC code or 
framework (none for parallelization, if the PIC code 
is already parallel). 



C. Free electron lasers 



In a short wavelength free-electron laser, a high 
energy electron beam interacts with a static mag- 
netic undulator. In the optimal boost frame with 
Lorentz factor 7, the down-shifted FEL radiation and 
up-shifted undulator have identical wavelengths and 
the number of required time-steps (presuming the 
Courant condition applies) decreases by a factor of 
27 2 for fully electromagnetic simulations. Exam ples 
of boosted-frame simulations have been compared [1 81 ] 
to results obtained with the eikonal (Le, SVEA) and 
wiggler-period averaged code Ginger [1 91 ] . It was con- 
cluded that if the necessary FEL physics can be stud- 
ied with an eikonal code, it will run much faster than 
a full electromagnetic code in whatever frame. How- 
ever, if there are important physical phenomena that 
cannot be resolved properly by an eikonal code, a 
boosted-frame electromagnetic code is a very attrac- 
tive alternative to a brute force full electromagnetic 
calculation in the laboratory frame. 



D. Coherent synchrotron radiation 

Another application for which the Lorentz-boosted 
frame method might be useful is that of modeling 
coherent synchrotron radiation (CSR) 20] emitted 
by high current, high brightness relativistic electron 
beams. Because full scale electromagnetic simulation 
of CSR in the laboratory frame is difficult due to the 
wide range of scales (chicane lengths of order meters, 
radiation wavelengths of orders microns), in order to 
make the calculation tractable most CSR simulation 
codes apply simplifications such as ignoring transverse 
variation of CSR across the electron beam. We have 
begun preliminary work of simulating CSR emission 
with the boosted frame method with Warp, examin- 
ing the behavior of a high current, short electron beam 
transiting a simple dipole magnet. Our early results 
show that upon exit from the undulator the electron 
beam shows the characteristic energy loss variation 
with longitudinal position that one expects from pre- 
vious theoretical analyses of CSR. Further studies are 
currently underway. 



IV. CONCLUSION 

The non-invariance of the range of scales of a phys- 
ical system implies that the computational cost of a 
certain class of computer simulations depends strongly 
on the choice of the simulation frame of reference. Al- 
gorithmic difficulties arise due to the loss of covariance 
upon discretization of the Maxwell- Vlasov system of 
equations, and the need to transform input/output 
data between the laboratory frame and the Lorentz 
boosted frame. So far, the difficulties that have arisen 
have been overcome and no "show-stopper" has been 
identified at this time. First principles simulations 
in boosted frames have been performed successfully 
with the code Warp in application to laser wake- 
field acceleration, electron cloud driven instabilities 
and free electron lasers, with speedups ranging be- 
tween a few and several orders of magnitude. Our re- 
cent progress show that first principles modeling in a 
Lorentz boosted frame is a viable alternative or com- 
plement to using reduced descriptions like the qua- 
sistatic [ll| or eikonal [l9| approximations, or per- 
forming simulations with scaled parameters [l7| . and 
in many cases includes physics that is not accessible 
to the other descriptions. This includes direct three- 
dimensional simulations of laser wakefield accelerator 
stages at 10 GeV and beyond, electron cloud effects 
in high energy physics accelerators, physics that is in- 
accessible to standard free electron lasers codes, and 
coherent synchrotron radiation. 



4 



Proceedings of the DPF-2009 Conference, Detroit, MI, July 27-31, 2009 



Acknowledgments 

We are thankful to D. L. Bruhwiler, J. R. Cary, 
B. Cowan, E. Esarey, M. A. Furman, C. Huang, S. F. 
Martins, W. B. Mori, B. A. Shadwick, C. B. Schroeder 
and M. Venturini for insightful discussions. 



This Work was supported by US-DOE Contracts 
DE-AC02-05CH11231 and DE-AC52-07NA27344, 
US-LHC program LARP, and US-DOE SciDAC pro- 
gram ComPASS. It Used resources of NERSC, sup- 
ported by US-DOE Contract DE-AC02-05CH11231. 



[1] J.-L. Vay, Phys. Rev. Lett. 98 (2007) 130405. 

[2] Proc. International Workshop on Electron- Cloud Ef- 
fects, Daegu, S. Korea (2007). 

[3] N. Kroll, P. Morton, M. Rosenbluth, IEEE J. Quan- 
tum Electron QE-17 (1981) 1436. 

[4] T. Tajima and J. M. Dawson, Phys. Rev. Lett. 43 
(1979) 267. 

[5] J.-L. Vay, Phys. Plasmas, 15 (2008) 056701. 

[6] J. P. Boris, Proc. Fourth Conf. Num. Sim. Plasmas, 
Naval Res. Lab., Wash., D. C. (1970) 3-67. 

[7] M. Karkkainen, E. Gjonaj, T. Lau, T. Wei- 
land, Proc. International Computational Accelerator 
Physics Conference, Chamonix, France (2006). 

[8] J.-L. Vay in preparation. 

[9] D. P. Grote, A. Friedman, J.-L. Vay, I. Haber, AIP 

Conf. Proc. 749 (2005) 55. 
[10] J.-L. Vay et al., Proc. Particle Accelerator Conference, 

Vancouver, Canada (2009) TU1PBI04. 
[11] C. Huang et ai, J. of Comput. Phys. 217 (2006) 658- 

679. 



[12] B. Feng, C. Huang, V. Decyk, W.B. Mori, P. Muggli, 
T. Katsouleas, J. Comput. Phys. 228 (2009) 5340. 

[13] C. G. R. Geddes et al, J. Phys. Conf. Series V 125 
(2008) 12002/1-11. 

[14] |http://loasis.lbl.gov| 

[15] D. L Bruhwiler et al., Proc. 13th Advanced Accelerator- 
Concepts Workshop, Santa Cruz, CA (2008) 29. 

[16] E. Cormier-Michel, et al. Proc. 13th Advanced Accel- 
erator Concepts Workshop, Santa Cruz, CA (2008) 
297. 

[17] C. G. R. Geddes et al., Proc. Particle Accelerator 
Conference, Vancouver, Canada (2009) WE6RFP075. 

[18] W. M. Fawley, J.-L. Vay, Proc. Particle Accelerator- 
Conference, Vancouver, Canada (2009) WE5RFP029. 

[19] W. M. Fawley, LBNL Tech. Rpt. LBNL-49625-Rev. 1 

(2004) ; also SLAC Rpt. LCLS-TN-04-3. 

[20] M. Venturini, et al, Phys. Rev. ST Accel. Beams 8 

(2005) 014202. 



